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We use lattice formulation of theory in order to investigate non-perturbative features of its 
continuum limit in two dimensions. In particular, by means of Monte Carlo calculations, we obtain 
the critical coupling constant g /in the continuum, where g is the unrenormalised coupling. Our 
final result is g/= 11.15 ± O.OGsiat ± 0.03si,st. 

PACS numbers: 12.38.Gc, ll.lSHa 


INTRODUCTION 

(j)^ theory plays a phenomenological role as an ex¬ 
tremely simplified model for the Higgs sector of the Stan¬ 
dard Model. In [HE] the triviality of 0^ theory in more 
than four dimensions has been proven, and there are nu¬ 
merous analytical and numerical results for D = 4 p[5], 
indicating that in this case the theory is trivial as well. 

In Zl = 2 and D = 3 the theory is super- 

renormalisable: the coupling constant has positive mass 
dimensions. In this paper we will work in D = 2, employ¬ 
ing lattice regularisation. In Z? = 2, [g] = [pq], where po 
is the (bare) mass parameter of the theory. This means 
that the only physically relevant dimensionless parameter 
is the ratio g/p^, where g is the bare coupling constant 
and p^ is a renormalised squared mass in some given 
renormalisation scheme. An additive mass renormalisa¬ 
tion is required since in the continuum limit the bare 
mass parameter diverges like log(a), where a is the lattice 
spacing. We do not care about coupling renormalisation, 
since it amounts to a finite factor. 

Despite the simplicity of the model, there is still de¬ 
bate in the literature about the value of / = gjp^, where 
the ratio is evaluated at the critical point. In particular 
we are interested in the value of /, call it /o, computed 
in the limit in which both g and p'^ go to zero; this corre¬ 
sponds to the critical value in the continuum. We decided 
to tackle this problem by using the same renormalisation 
scheme used in lilll, adopting the simulation technique 
introduced in |5], namely the worm algorithm, and us¬ 
ing a completely different strategy to obtain g/p^ in the 
infinite volume limit. 

In the following we will describe the model and the 
renormalisation scheme chosen in order to extract p^ at 
fixed g in the infinite volume limit from our simulations. 
Then we will give details about the simulations and 
we will proceed to the continuum limit extrapolation. 
In the end we will compare our results with recent 
determinations of the same quantity and we will draw 
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some conclusions. 


I. LATTICE FORMULATION 

Let’s introduce the (j)'^ Lagrangian in the Euclidean 
space: 

= \ -b (1) 

In ZZ = 2 the Euclidean action is 



In order to obtain a dimensionless discretized action we 
put the system on a 2-dimensional lattice with spacing a 
and introduce the following parametrization 

Fo = a Pq, g = a g. (2) 

In this way we have 

= E 1“ E + 2 (^0 + 4) (fi + : (3) 

where ((>x±o are fields at neighbor sites in the Ini' direc¬ 
tions. 

In the following we will omit the “hat” on top of lat¬ 
tice parameters: all quantities will be expressed in lattice 
units, i.e. they become dimensionful when multiplied by 
appropriate powers of the lattice spacing a. 

If we take the continuum limit too naively, at fixed 
physical quantities, we obtain, in ZZ < 4, the critical 
Gaussian model (5). On the other hand, if we stick to 
a fixed value of g (in lattice units) we can search for a 
value of Pq such that we get, in the infinite volume limit, 
a second order phase transition point in the plane {g, Pq). 

In order to safely go to the continuum limit, we have 
to work out an additive renormalisation of the mass pa¬ 
rameter, since p^ in this limit diverges like log(a); in 
this way we translate Pq into p'^, a renormalised squared 
mass. Of course several definitions of renormalised mass 
can be chosen; in this work we adhere to the same renor¬ 
malisation procedure as in 13 El- We refer the reader to 
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these papers for more details. Here we only remind that 
in D = 2 there is only a 1-Particle-Irreducible divergent 
diagram (see Fig. [^. Its expression on a lattice with 
N X N points is 

Af-l N-l 

ki—0 k2—0 4 

( 4 ) 

and a suitable renormalisation condition consists in 
putting equal to the solution, in the infinite volume 
limit, of the equation 

=^il + 3gA{^i^). (5) 


sm 


Trfci 

w 


sm 


T:k2 

Iv 


Mo 



This condition is equivalent to the introduction of a 
proper divergent mass-squared counterterm in the ac¬ 
tion. We may finally extrapolate the quantity / = 
to g —^ 0 in order to obtain /o, the critical value in the 
continuum limit. 

Another parametrization of the action is the following: 

^ ^ + ^Wi + - 1 )^] 

X . X (6) 

= Si + Ssite-, 

where the relations between (/Tq, g) and (/3, A) are: 

<(>x = nl=2— -4, 9=^- ( 7 ) 

In eq.([^ there is an interaction term between neigh¬ 
bor sites, Si, with a coupling constant of strength /3 
and a term related to a single site, Ssue- With this 
parametrization it is easy to recognize the Ising limit for 
A —^ 00 . In this limit, configurations with ^ 1 are 
completely suppressed and the fields assume only values 
'p{x) = ±1. As a result, the second term of can be 
disregarded and the action becomes the well-known Ising 
action Se = -/3 Y.u Vx>Px+o- 


A. Simulations 

In this section we outline our general computational 
strategy, postponing the discussion of the simulations de¬ 
tails. 

We use the worm algorithm [8] , using the lattice action 
given by Q. We checked our simulation program against 
the results of mm, obtaining values compatible within 


errors, well below one sigma level. In this case and also 
in the following, in order to estimate statistical errors we 
use the program described in El- 

Considering a fixed value of A, our aim is to compute 
the critical point of the theory, i.e. the critical value of 
/3 for that particular value of A. We use the physical 
condition 


mL = L/^ = const = 2 , (8) 


where m is implicitly defined by the condition 

G{p*) wS' 

G'(O) -|- 


(9) 


G{p) is the two-point function in momentum space, and 
p* is the smallest possible momentum on a lattice of lin¬ 
ear size L. Details, as before, in [8]. Condition Q im¬ 
plies that ^ grows linearly with L, and when L/a —>■ oo 
we arrive at the critical point. We then simulate several 
lattices with different values oi N = L/a; for each couple 
(A, N) we obtain a value of /3(A, N) such that mL = z. 
After this step we extrapolate our results to a/L —>■ 0 in 
order to compute /3(A). Now, using relations in 0 we 
derive g{X,l3) and Pq{X,(3). Using renormalisation con¬ 
dition we finally pin down g^{g) and hence the ratio 

/ = g/h^- 

We repeat all this procedure for several values of A, 
and hence of p; in the end we extrapolate our results to 
p —>■ 0, in order to obtain /q. We will now focus on the 
details of our simulations. 


We choose the condition z = 4. As we will see in the 
following, this choice is not as crucial as it may seem. 

At a fixed value of A we simulate the system for five 
values of L/a, namely: L/a = 192, 256, 384, 512 and 
768. For each value of L/a few preliminary simulations 
are needed to roughly find the value of /3 leading to z ~ 4. 
In few cases (see for example Fig. we have explicitly 
checked that using five values of /3 such that z falls ap¬ 
proximately into the interval [3.8,4.2] we do not observe 
any sign of non-linearity of z as a function of j3. The 
difference in ,0(z = 4) between the case in which we use 5 
points to interpolate and the case in which we use only 3 
points is one order of magnitude less than the statistical 
error itself. We then decided to use just 3 values of /3 
for the real simulations to linearly interpolate the results 
and to obtain in this way /3(A, A). 

A typical full simulation (A = 0.25) is synthesized in 
Table |T] is the number of worm-sweeps between 

two measures, which increases in order to minimize the 
simulation time, taking into account autocorrelation 
time; the number of thermalisation sweeps for all our 
simulations is several hundreds times r, the autocor¬ 
relation time of mL, which we always keep under control. 

(j)'^ theory [12] is in the same universality class of the 
Ising model, and we know that in D = 2 the critical expo¬ 
nent of the correlation length is ^ = 1. Thanks to finite 
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2. Linear interpolation (A = 0.25, L = 

256) to obtain 

d(z = 
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TABLE I 

A = 0.25 simulations 


L/a 

^mens 

sweep 

fidz = 4) 

192 

1 X 

15 

0.655357(12) 

256 

5 X 10"^ 

15 

0.656177(11) 

384 

5 X lO'^ 

15 

0.656984(8) 

512 

3 X 10"^ 

20 

0.657399(7) 

768 

2 X 10"^ 

25 

0.657818(10) 

OO 



0.658628(10) 


dent on the reliability of the extrapolations. As an ex¬ 
ample we show, in Fig. a double extrapolation to 
a/L = 0 in the case A = 1. For z = 4 the extrapolation 
to a/L = 0 is steeper than for z = 1, since in the latter 
case, at finite volume, we are nearer to criticality, so that 
/3(A, N) is not so far from the infinite volume value. Nev¬ 
ertheless at z = 4 we obtain a much more clear signal; 
we can extrapolate to the a/L = 0 value with a much 
smaller statistical error even if the number of measures 
is (5 — 10)-times smaller than the case z = 1. The results 
in the infinite volume limit coincide within the statistical 
errors; /3(z = 1) = 0.68060(4), to be compared with the 
equivalent value in Table [III = 4) = 0.680601(11). 
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0.680 

0.679 
^ 0.678 
0.677 

0.676 

0.675 

0.674 
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size scaling arguments we expect to be able to extrapo¬ 
late P{X,N) to (3{X) linearly in a/L. This is numerically 
very well confirmed for all values of A we explored. In Fig. 
[^we show a typical extrapolation. For every value of A 
considered, we obtain a very reasonable value of ^ 1- 
Our final results are reported in Table |TTj 


FIG. 4. Extrapolation to a/L — 0 with z = 4 (blue steep 
curve) and z = 1 (red curve) (A = 1). 


II. RESULTS 
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FIG. 3. Linear extrapolation of /3 to a/L — 0 for A = 0.25. 



Now we show that the condition z = 4 is not cru¬ 
cial; actually, as is well known from general theoretical 
arguments, we could choose another value of z without 
affecting the results in the infinite volume limit. From 
a numerical point of view it is nevertheless interesting 
to consider other values of z in order to be more confi- 


In Fig. [^we plot the results shown in Table [Il| The plot 
is in x-log scale, to emphasize the fact that we covered 
over two order of magnitude in g. Blue round points are 
our results taken from Table [n| Red triangular points are 
results from [7] . We postpone the discussion of the green 
square points. 

First of all we note that in the intermediate region, i.e. 
in the minimum of the curve, our results are in almost 
perfect agreement with those of [7] . Note that the infinite 
volume limit results of [7] are obtained with a completely 
different strategy. The situation starts changing at the 
lowest simulated values oi g: we see, in the insert shown 
in Fig. that our points seem to be a little bit higher. 
The blue curve is our final fitting function, which we are 
now going to discuss, while the red dashed curve is the 
fit function used in [7]. 

We decided to fit f{g) over the entire range at our 
disposal with the function 


fid) 


gp -I- aig + a2g^ + + Q4g^ 

1 + 


( 10 ) 


We can certainly justify the functional form for large val¬ 
ues of g. We know that (j)'^ theory reduces to the Ising 
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TABLE II. Final extrapolations to infinite volume limit: g and g? are computed at Pc using equations ^ and 0. 


A 

Pc 

g 

h 

gig 

1.000000 

0.680601(11) 

8.63523(29) 

0.649451(67) 

13.2962(18) 

0.750000 

0.689117(13) 

6.31733(24) 

0.509730(59) 

12.3935(19) 

0.500000 

0.686938(10) 

4.23833(12) 

0.367173(31) 

11.5431(13) 

0.380000 

0.678405(11) 

3.30267(10) 

0.296195(32) 

11.1503(15) 

0.250000 

0.6586276(98) 

2.305261(69) 

0.214762(27) 

10.7340(17) 

0.200000 

0.6462478(78) 

1.915543(46) 

0.181077(21) 

10.5786(15) 

0.125000 

0.6190716(52) 

1.304633(25) 

0.125924(15) 

10.3605(15) 

0.094000 

0.6030936(89) 

1.033757(30) 

0.100518(23) 

10.2843(26) 

0.062500 

0.5820989(60) 

0.737813(15) 

0.072073(15) 

10.2370(23) 

0.030000 

0.5516594(71) 

0.394311(10) 

0.038407(17) 

10.2666(48) 

0.015625 

0.5326936(27) 

0.2202547(22) 

0.0211916(63) 

10.3935(32) 

0.007500 

0.5187729(29) 

0.1114722(12) 

0.0105457(67) 

10.5704(68) 

0.005000 

0.5136251(17) 

0.07581192(49) 

0.0071014(38) 

10.6757(57) 

0.002000 

0.5064230(16) 

0.03119343(19) 

0.0028637(35) 

10.8925(132) 



FIG. 5. Final results for f{g) in logarithmic scale. Error-bars, 
where not visible, are smaller than symbols size. 


model in the limit A —> oo. In particular in the Ising limit 

we have /? = ^ . Note that /3(A), at 

the critical point, is a highly non-linear function of A 
itself. In fact at A = 0, /3 = 0.5; then we note a maxi¬ 
mum, with a value around 0.69 for intermediate values of 
A; in the end /3(A) has to go asymptotically to the value 
0.44068679 ..., the critical Ising value in H = 2. In [T3] 
it is noted that for A = 10 the value of /3 at criticality is 
already near the asymptotic value. For very large values 
of A we can then safely approximate P with if we 

look at the relations Q, we note that g is going to in¬ 
finite linearly with A, and /ip diverges proportionally to 
g. But this is not true for g? due to the renormalisation 
condition (©• We numerically checked that gP, using the 
approximation /3 = for g > 10^, can be linearly ex¬ 

trapolated in 1/g to g ^ oo (see Fig. [^. We arrive at 
the value g'jging = 3.40669(1); the error is subjectively 


estimated from the fit. 

We simply assume a linear behavior of f{g) for 5 —>■ 0. 
Taking into account the Ising limit constraint, we fix the 
parameter 63 as a constant times 04 . We have in total 7 
d.o.f. and we obtain 

fo = 11.179(62) (11) 

with a reduced = 0.73. 



FIG. 6. Extrapolation of g^ at g ^ 00 , as described in text. 


In order to check the validity of the ht function (10 1 , we 
decided to compute f{g) with the same strategy adopted 
in [7] , but for two values of g higher than those considered 
in [7], namely 5 = 4 and 5 = 6 . The field configurations 
are generated with a mixture of Metropolis steps and 
single cluster Wolff steps, used in [7] and presented in 

[la¬ 
in particular for each L/a we search for the value of 
/Tq that maximize the magnetic susceptibility x = w) ~ 
) ^; this peak is a signal of the pseudo-transition point 
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at finite volume, (p is the average of the held over the 
whole lattice. is then extrapolated to a/L —>■ 0 and 
the corresponding is obtained by means of condition 


(§. 


Details of simulations for g = 4 are given in Table III 


As can be seen in Fig. the two points at (/ = 4 and 


TABLE III. (? = 4 simulations with Metropolis-cluster algo¬ 
rithm 


L/a 

^meas 

glg^ 

128 

1 X lO'" 

11.2631(13) 

192 

1 X 10 ® 

11.3227(9) 

256 

1 X 10 ® 

11.3533(7) 

384 

1 X 10 ® 

11.3826(3) 

512 

1 X 10 ® 

11.3969(3) 

00 


11.4417(5) 


g = Q, represented by squares, lie perfectly on the curve 
dehned by our ht function. This represent a further con- 
hrmation that our strategy for computing g /passing 
through the limiting procedure described above, works 
as expected. 

In order to better understand the behavior of f{g) for all 
possible values of g we dehne a new parameter, 77 : 


g = 


5+1 


( 12 ) 


It is clear that (12) is a map from g S [0, 00 ) to g G [0,1]. 
We hope in this way to obtain a smoother behavior of 
f{g); note that the limit f{g -G 0 ) is completely equiva¬ 
lent to f{g —)■ 0). We then dehne the ht function 


^ a'o + a[g + a' 2 g^ + a'^g^ 

1 + b[g + b'^g^ + b'^g^ ’ ^ ^ 

where one of the parameters is determined by the Ising 
constraint for 77 = 1 . 

As shown in Fig. this choice leads us to a smoother 
function. With the g parametrization we obtain: 


/o = 11.119(24), 


with a reduced = 0.95 and 8 d.o.f. 


(14) 


III. CONCLUSIONS 


We decide to quote our hnal result as: 

/o = 11.15(6)(3). 


(15) 


We take as central value the mean of 0 and ( [T4| ). The 
hrst error is purely statistical, and it is conservatively 
taken as the biggest one between the two hts. The second 
error is an estimate of the systematic error associated 
with the particular functional form used to ht data. 

In Table ITVl we summarize some of the latest results for 
/o derived with different approaches: the works [HHni 



FIG. 7. Final plot of f{g) with our results. 


Uni |20] are based on Hamiltonian truncation (variational) 
methods, while in [TH] lattice theory is simulated by using 
non-local SLAG derivative. 

We note that our result is compatible with the last 
four determinations, which come from different methods. 
We only observe a discrepancy at a Scr-level with the 
Monte Carlo results in [7], where a region of very small 
( 7 -values is reached. For technical reasons, which will 
be hopefully overcome in the near future, we could not 
reach this region, but thanks to the worm algorithm our 
statistical errors are much smaller. We also note that the 
result of our second fit ( 77 -parametrization, see Fig. 
has a statistical error comparable with that of m, and 
the two results are compatible at 2CT-level. Although 
we were very conservative in the error estimations, we 
believe that this work is a step towards a more precise 
Monte Carlo determination of /q. 

Our plans for the next future are to improve this work 
towards the g —>■ 0 limit with an extended statistics. 


TABLE IV. Sample of the results for the continuum crit¬ 
ical parameter /o from the literature. DLCQ stands for 
Discretized Light Cone Quantization, QSE diagonalization 
for Quasi-Sparse Eigenvector diagonalization and DMRG for 
Density Matrix Renormalization Group 


Method 

fo 

year, Ref. 

DLGQ 

5.52 

1988, [IS] 

QSE diagonalization 

10 

2000, [IS] 

DMRG 

9.9816(16) 

2004, [H] 

Monte Garlo cluster 

10.8g:),5 

2009, 13 

Monte Garlo SLAG derivative 

10.92(13) 

2012, [IS] 

Uniform Matrix product states 

11.064(20) 

2013, [E] 

Renormalised Hamiltonian 

11.88(56) 

2015, [IS] 

Monte Garlo worm 

11.15(6)(3) 

This work 
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